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DESIGN OF A HELICOPTER AUTOPILOT BY MEANS OF LINEARIZING TRANSFORMATIONS 
G. Meyer, R. L. Hunt, and R. Su 

NASA Ames Research Center, Moffett Field, California, U.S.A., 94035 


SUMMARY 

A method for designing automatic flight control , systems for aircraft that have complex characteristics 
and operational requirements, such as the powered-lift STOL and V/STOL configurations, is presented. The 
method is effective for a large class of dynamic systems that require multiaxis control and that have highly 
coupled nonlinearities, redundant controls, and complex multidimensional operational envelopes. The method 
exploits the possibility of linearizing the system over its operational envelope by transforming the state 
and control. The linear canonical forms used in the design are described, and necessary and sufficient condi- 
tions for linearizability are stated. The control logic has the structure of an exact model follower with 
linear decoupled model dynamics and possibly nonlinear plant dynamics. The design method is illustrated with 
an application to a helicopter autopilot design. 


1. INTRODUCTION 

Consider in general terms the control -system design problem. Let us take the usual hardware/software 
model of the problem in which the hardware consists of the plant together with all the sensors and actuators 
which are connected to a digital computer, and in which the software embodies the complete control strategy. 
The hardware is fixed; we may change only the software, over which, however, we have full control. So one 
may say that since the underlying physical process is to remain fixed, only its representation may be changed. 
If this point of view is taken, then much of the control-system design problem may be interpreted in terms 
of transformations. 

This paper outlines a design approach that is being developed from the transformations point of view, 
and describes an application to the control of a helicopter. This approach, first outlined in Ref. 1, has 
been applied to several aircraft of increasing complexity, and the completely automatic flight-control system 
was first tested on a DHC-6. The reference trajectory used in the flight test .exercised a substantial part 
of the operational envelope of the aircraft. Despite disturbances and variations in plant dynamics, the 
system performed well (see Ref. 2). Next, the technique was applied to the Augmentor Wing Jet STOL Research 
aircraft, the successful flight tests of which are reported in Ref. 3. Methods for providing pilot inputs to 
this design were examined in Ref. 4, and application of the scheme to the controV of an A-7 aircraft for 
carrier landing and testing in manned simulation is reported in Refs. 5 and 6. The design method is currently 
being applied to the UH-IH helicopter, again with the substantial portion of the operational envelope of this 
aircraft being used. 

The key concept of the approach is to simplify the representation of the plant dynamics by means of a 
change of coordinates of the state and control. The design proceeds in three steps. First, the given non- 
linear system - possibly time-varying, multiaxis, and cross-coupled - is transformed into a constant, decou- 
pled linear representation. Second, standard linear and nonlinear design techniques, such as Bode plots, pole 
placement. LQG, and phase plane, are used to design a control law for this simple representation. And third, 
the resulting control law is transformed back out into the original coordinates to obtain the control- law in 
terms of the available controls. 

The mathematical foundation for the approach is provided by modern differential geometry, the necessary 
and sufficient conditions for linearizing have been established, and the theory is given in Refs. 7-15. 


2. CANONICAL FORM 

Suppose that the given natural representation Si of the physical process is given by the state Xi, 
control Ui, and field fi, 

Xi = fi{Xi,Ui) , (1) 

and that we wish to change Si into S 2 with state X 2 » control Uz* snd field fz, 

X 2 = f2(X2,U2) (2) 

This change will be accomplished, and hence a large part of the design problem solved, once we construct the 
appropriate transformations of the state and control, 

xz = T(xi) (3a) 

= W(xi,U2) (3b) 

which relate Si to S 2 so that for all admissible (xi.Uz), 

(xx)fi[xi,W(xi,U:)] = f2[T(xx),U2] ( 4 ) 

The function W is the control law. 

The construction of the transformation (T,W) is often greatly simplified by the introduction of an 
intermediate canonical representation Sq ss shown in Fig. 1. To obtain (Ti 2 *Wi 2 )j which links Si to 02 * 
both Si and S 2 are first transformed into the canonical representation Sq* Then 
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Ti 2 - TiiTio 

Wi2 - WioWJ} 



Fig. 1. Manipulation of system representations. 


»o m • • ek^ INTEGRATORS 

^2 >• • • • • • Kj INTEGRATORS 

”o ,,,,, 

uy • • . • • INTEGRATORS 

Fig. 2. Diagram of Brunovsky form. 


(5a) 

(5b) 


In the design procedure being described, the 
Brunovsky form (Ref, 14) Is taken to be canonical. 

In linear theory this form Is basic. It consists of 
decoupled strings of Integrators which may be dia- 
granmed as shown In Fig. 2 where each dot represents 
a scalar Integrator. 

The number of strings, which may be of different 
lengths, equals the number of controls which Itself 
equals the number of Kronecker Indexes 
ki i k 2 a . . . > kin. The lengths of the 1th string 
Is given by ki and the dimension of the state space 

m 

n = 52 ^4 • 

1=x ’ 


Let the canonical state Xq e r” con .rol Uq e r'”, 
and denote the canonical field fo by 


0 0 0 0 


( 6 ) 


According to linear theory any constant, linear, 
controllable system may be viewed as a nonsingular 
transformation of an appropriate Brunovsky form. 

Thus, if Si is given by the state Xi g R", control 
Ui G R^, and field 


Xi * AiXi + BiUi 


(7) 


with controllable (Ai,Bi), then there Is a Brunovsky form that can be transformed Into Si by means of non 
singular transformations (T,W nonsingular) 


Xl ■= T-‘Xg 


(8a) 


Ui * WUq + Qxi 


(8b) 


Of course, S© nay also be transformed Into nonlinear and time-varying systems without any loss of 
Information by allowing the transformations In Eq. (8) to be nonlinear and time-varying but still nonsingular: 

Xl = T’Mxo.t) (9a) 

ui » W(xi,UQ,t) (9b) 


This fact gives rise to our design procedure. When presented with a nonlinear system Si, the first step 
Is to try to linearize the system over Its whole operational envelope by constructing (T,W) which maps Si 
Into Sq. Then a control law Is synthesized for the much simpler S©. Finally, the S© control law Is 
transformed Into the coordinates of Si to obtain the control law for Si. 


3. TRANSFORMABILITY 



Fig. 3. An example of a block triangular system. 


A class of systems particularly amenable to this 
approach has the following form. The control u 6 Rf^, 
the state x g Rf" x R'" »< . . . x R*” and the field f 
is without transmission zeroes and invertible. An 
example is shown in Fig. 3. In this example, the 
state x G R"’ X R"’, the control u G R"’, and the 
field f. 


Xl - Fi(xi,X2,t) 

X 2 » Fi(xi,X 2 *u,t) 


( 10 ) 


and Fi and F 2 are Invertible relative to (x 2 ,xi) and (u,X 2 ). That Is, functions hi and hz can be con- 
structed so that If 

X2 - hi(Xi,Xi.t) 

U « h-(Xi,X2,X2,t) 



then 


Fi[xi,hi(xi,xi,t),t] - Xl 

F 2 CX 1 »X2 »h2 (Xi ,X2 »X2 *t) ,t] “ X 2 


( 12 ) 
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on the operational envelope. Because of the form of Eq. (10), such systems will be called block- 
triangular. 

For this example, the canonical form Sq has m Kronecker indexes, all equal to 2. The state 
xo^ X Rtri, the control u° e R^, and the field, 

X? = x! j 

xl = u'’ j 


(13) 


The map linking Eq. (10) with Eq. (13) may be obtained by letting Xi = x?(t) and pushing, as it were, 
the time-history Xi(t) upstream through f to obtain u(t). Thus, 

xi = x! \ 

xi = X? = x! I 


X 2 = hi(x?.x!.t) 


X2 


_ 8hi x° + 


9Xi 


9X2 


9t 


(14) 


U = h2(Xi ,X?,X2,t) 


In general, the canonical form Sg of a block triangular Si will have m Kronecker indexes, all 
equal to n/m where n is the dimension of the state space of Si. 

This procedure of constructing the linearizing transformations will fail if Si does have transmission 
zeroes. In that case, one obtains differential equation constraints on u, thereby destroying its status as 
an independent control variable. Nevertheless, such systems may still be lineanzable. For example, the 
scalar system 



(15) 


is linearized by the transformation 

xf - Xi + i (xs)^ 


Xf = X2 
X? = X3 
u = 


(16) 


On the other hand it will be shown that the system 



is not linearizable; the nonlinearity in this case is intrinsic and cannot be removed by a change of coordi- 
nates. The conditions under which linearization is possible are summarized next. 

Let X e u 6 R^, and the field 

X = F(x,u) (1®) 

There are four conditions for this system to be linearizable. First, it is necessary to be able to construct 
a new control variable, 


= h’^(x,u) 
u = h(x, 4 >) 


(19) 


so that <{) enters linearly into the field: 

m 

F[x,h(x,(())] = f(x) + S gi(x)<i>j 
i=i 


(20) 


The remaining three conditions are technical, and they are best expressed by means of Lie brackets defined as 
follows. If f and g are C” vector fields on R"^, the Lie bracket of f and g is 


[f.g] = f f - 


9X 


( 21 ) 



4 


and we set 

(ad“f,g) = g 
(ad*f,g) = [f,g] 
(ad^f.g) = [f,[f,g]] 


(ad'^f.g) • [f.(ad'^-‘f,g)] 

A collection of C” vector fields hi.hj, . ...hr on R" is involutive If there exist C” functions 

^ijk 

r 

[h^.h.] » 2 Y^j|jh,j , Isi, j<r, 1)<j 
k*i 

Now, suppose that we wish to transform Eq. (20) into Sq with Kronecker indexes ki > k: > . . . > km. 
Define the sets 

C = {gi.[figi]. • . •» (ad*^^ ^figi)* 
g2»[^ig2] (ad^^ ^fig2)i 

v[7xiL'“.“.T7a?"*"f7i)) 

Cj = (gi.Cfigi]. .... (ad*''5*^figi). 

gj.[fi9i]» • • •• (ad^j fig2). 

9ra*t^‘9ml* ' * (ad‘‘7''fig^„)) for j = 1 . 2. . . .. m 

Then it can be shown that the transformation is possible if and only if at each admissible x, 

1. The set of C spans an n-dimensional space 

2. Each Cj is involutive for j = 1, 2, . . m 

3. The span of Cj equals the span of Cj n c for j = 1, 2, . . m 

For a linear field, x = Ax + Bu, the spanning condition (1) or C is equivalent to controllability, rank 
(B, AB, . . ., A>^"^B) « n. The other conditions are automatically satisfied. The new coordinate surface 
Ti(x) = constant in a plane through the origin in the old state space. For nonlinear field. Ti(x) = constant 
will be a general surface. The involutivity condition (2) guarantees that this surface is constructible from 
local conditions (integrability theorem of Frobenius). 

Suppose we wish to transform system (17) into Sq with k = 3. Here the set C = {g,[f,g],(ad^f,g)} 
spans R^ But the Ci * {g,[f ,g],[g.[f .g]]> also spans. therefore, system (17) is not transformable into 

Sq. For further details on the transformation theory see Refs. 7-15. Let us turn our attention next to the 
control structure in which these ideas of transformability may be implemented in practical cases. 


4. EXACT MODEL FOLLOWER 

If the plant $i is equivalent to a system S 2 in the sense that each can be transformed into the other 
by nonsingular transformations, then one may construct an exact model follower in which the plant Si will, 
except for disturbances, follow exactly the system $ 2 * which is Interpreted as being the model. Let the 
system representing the tracking error be denoted by S3. The three systems, are related to each other and 
to the Brunovsky form Sq as shown in Fig. 4. 



Fig. 4. Plant Si, model S 2 , and regulator S3. 
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Now consider the structure of the model fol- 
lower as shown in Fig. 5. There are five subsys- 
tems; model servo, regulator, plant, and two 
transformations.’ The desired system behavior is 
defined in the model servo. The model servo law 
may be nonlinear, time-varying, and dynamic. 

The TW-map (T 23 ,W 23 ) transforms the S 2 system 

xT = f2(x?,u?) . (22) 

into S3. 

X? = fj(x'?,u'!I) (23) 

The TW-map (Tis.Wu) transforms the plant Si 



Xi = fi(xi,ui) (24) 


Fig. 5. Structure of the model follower. 


into S 3 


X 3 - f3(X3.U3) (25) 

That is, when seen through the TW-map (Tis.Wia) in termsof (X 3 ,U 3 ), the plant looks like Eq. (25); which is 
also how the model looks through the TW-map (T23.W23). 

In the absence of disturbances and with proper initialization of the model (i.e., X 2 ( 0 ) = Ti 2 [xi( 0 )]), 
the plant will follow the model exactly in the sense that the error defined by . . : 

63 = X 3 - X? ' (26) 

will be zero for t > 0. In a realistic situation, disturbances, say, d, are present, and they are controlled 
by means of the regulator which transforms the tracking error 63 into corrective control 6 U 3 . The error 
dynamics are given by 

63 = f3(X3»U3) - T3 (x?»U3*) (27) 

which for small errors simplifies to 

In particular, if S3 is chosen to coincide with Sq. so that fs = fo, then without approximation, 

^0 = AgBo + B(j6Uo + d , (28) 

It may be noted that the regulator need not be gain-scheduled; that function is accomplished automatically 
by the TW-map. 

This completes the outline of the design approach. Consider next an application to a helicopter. 


5. THE PLANT-A HELICOPTER 

The helicopter will be represented by a rigid body moving in three-dimensional space in response to 
gravity, aerodynamics, and propulsion. The state, 

X = (r.v.C.u)^ e X C R’ X R5 X S0(3) X (29) 

where r and v are the inertial coordinates of body center-of-mass position and velocity, respectively, and 
C is the direction cosine matrix of the body-fixed $xes relative to the runway-fixed axes (taken to be 
inertial). The attitude C moves on the sphere S0(3). The body coordinates of angular velocity are 
represented by (jj. . , . , . , 


The controls, 

u = (u^.u’’)^ e uc R= X R (30) 

where is the three-axis moment control, that is, roll cyclic and pitch cyclic, which tilt the main-rotor 
thrust, and the tail-rotor collective, which controls the yaw moment; and u^ is the. main-rotor collective, 
which controls the main-rotor thrust. 


The effectively 12 -dimensional state equation consists of the translational and rotational kinematic and 
dynamic equations: 


r = V 


V = f'^(x,u) 
t » S(u)C 
U=‘f”(x,u) /■ 


(31) 
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where and are the total force and moment generation processes, and 
and (30). This Is the natural representation Si of the helicopter. 


(x,u) are defined by Eqs. (29) 


Consider next the transformation of Eqs. (31) Into a Brunovsky form So by means of the transformations 
u = W(x,Uo) and Xq ■ T(x). 


6, THE TW-MAP 

In general, the moment generation process f^ Is Invertible with respect to the pair (u,u^), and for 
the restricted class of maneuvers being considered in this experiment (I.e., no 360® rolls), f^ Is 
Invertible with respect to the pair (vi,u*^). Thus, in the present case for the set of angular and vertical 
acceleration commands restricted to the set 

• {{i,V|) : 5 1.0 rad/secS 1 • 1,2,3, Ivj] $ 0.5 g) (32) 

a function h^:X U can be constructed so that If 

u « h”[x.(wo,v,Q)] 


then 


U « (JQ I 

S - ^30 I 


(33) 


(34) 


If (wq.Vjo) are chosen to be the new Independent control variables to replace the natural controls (u^.u^) 
then the state equation (31) becomes the following: 


■ f‘’(r,v,C) + c f^[r,v,C,o),((Lo,v,o)] 


(35) 


*30 

t « S(u))C 
u - Oq 

where e * 1 and f^ Is such that f'[r,v,C,0,(0,0)] » 0. 

The function f® Is Invertible with respect to the pair, {[vi,; 2 .E 3 (il»)],C} In which is an elemen- 

tary rotation about the runway z-axis, representing the heading of the helicopter. 

/ cos sin 0 \ 

Ej(\|>) » I -sin ri> cos 0 I (36) 

\ 0 0 1 / 

If the horizontal acceleration commands are restricted to the set 

■ ((vi.Vj) : |v^| s 0.5 g, 1 « 1,2} (37) 

then a function h^ : R* « R* * k S0(2) -► S0(3) can be constructed so that the helicopter attitude given by 

Cq « h^[r,v,VQ,EJ(,^»o)] (38) 


results In the comnanded acceleration. 


V ■ Vo (39) 

Equations (33) and (38) are the trim equations of the helicopter (31) without the parasitic effects (e = 0 
In Eq. (35)). That is, for a given motion {r(t),E,[ijr{t)]},t z 0. the corresponding trim state and control 
may be computed as follows: 

•■o “ Kt) 

Vo • r{t) 

Co “ h^rj,,Vg,Vj,(t).£,[»o(t)3) 

"0 - 0[Co(OcJ] 

«0 * “ot^) 

“o • h”t»-o*''o*Co*“o*<“o'''so)3 


( 40 ) 
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where the function q extracts u from S(tj) - Cc .. .The, 
required time derivatives in Eq. (40) are computable pro- 
vided that the motion (r.Ea) is generated by the strings 
of integrators shown in Fig. 6 where a dot represents a 
scalar integrator, and e R** is an independent control . 
variable. The system shown in Fig. 6 (Kronecker indexes 
{ 4 , 4 , 2 , 2 }) will be taken as the canonical model of the 
helicopter. The canonical variables will be denoted by'’ y 
rather than Xq to reduce the number of subscripts. The 
transformation and feedback that change the natural repre- 
sentation (Eq. (31)) to the canonical representation are 
approximately the following. 

The coordinate change, y = T(x), is given by 

y' = (y'l.yp^ ' 
y’ ' (yl.yl.yj.yj)^ = 


y" = (yJ.y^.yJ.yJ)"'' 


af" . T 
IT 


(41) 


where 


y 

•- 


Y 




-♦ Eo 






. Fig. 6 . Canonical model of helicopter. 


;COS i> = c^^/c, 'Sin ^ ~ ^ 

The control variable change, u = w{x,y®) is defined in two steps: 



(42) 


and 


u = h^[r,v,C,u,(o),V 3 )] 


(43) 


The effects of the various approximations made in the construction of these transformations are relegated to 
the regulator. 


7. MODEL SERVO AND REGULATOR . ' . ; 

The field of the model $2 is chosen to be canonical as defined in Fig.- 6 . The design of the three 
trajectory channels of the model servo is shown in Fig. 7. There are four 3 -axis integrators. Two axes 
correspond to the horizontal strings (ri,r 2 ) in Fig. 6 . The third axis Unla?* 

but with two additional smoothing integrators. The exceedingly simple structure of the field (i linear, 
decoupled) greatly simplifies the servo design 'process. The inner loop (ao,ao) is designed to be an acceler- 
ation servo^whose input ai ’ is the sum of coarse acceleration command a* from the coarse command gener- 
ator, and 6 ao wTiich is generated by the outer loop to smoothly reduce any 

position and velocity vectors, r* and v*, respectively. The acceleration servo bandwidth is 0.63 ^^d/sec. 

The bandwidth of the outer loop is 0.1 rad/sec. Large position errors are reduced at the rate of 6 m/sec. 

The heading model servo is shown in Fig. 8 . It is designed to have two scalar integrators corresponding 


to the Ej .string in Fig. 6 . The heading error, 


cos tpQ sin i)>* - sin cos ip* 


(44) 


is computed in the q-block. Block S represents, together with the integrator, the kinematic equation 
of E 3 . 

The model reference state yo and control y^ in Fig. 5 are defined (see Figs. 6 - 8 ) as follows: 


yj = (»'io*’^2o) 
yj ” ('^io»'^2o^ 

y© * fa^Q »ajQ jEj ( i^q) ,rgQ] 

y© " (^l0*®20»‘*^30»^30^ 
yo ■ (^*10»^20*‘^30»®30)*^ 


T . 


(45) 



INITIAL 

WAYPOINT 


r(o) 


COARSE v*+t- 

TRAJECTORY^ Mi— 

GENERATOR 


v(o)— » j / 
^ I 

S T 






Fig. 7. Model servo-trajectory, r « (rj,r2,rj). 



Fig. 8. Model servo-heading £3(11^0). 


The field of the regulator S, is also chosen to be canonical Sq. The regulator design is outlined 
in Figs. 9-11. 


Fig. 9. Regulator- horizontal axes. 


Fig. 10. Regulator - vertical channel 
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The estimated canonical state y is computed 
from the estimated natural state by means of the 
transformation T defined by Eq. (41). The outer 
loop of the horizontal channels (Fig. 9) has a 
bandwidth of 0.3 rad/sec, with a 3 m/sec large- 
error reduction rate and 0.1-g authority limit. 

The inner loop bandwidth is 3 rad/sec. The verti- 
cal channel, shown in Fig. 10, has a bandwidth of 
0.63 rad/sec, a 1.5-m/sec large-error reduction 
speed, and 0.3-g authority. Finally, the heading 
channel is regulated as shown in Fig. 11, where 
the heading error is given by 

e. = cos 4; sin 4;- - sin 4> cos 4/ (46) 

ij; ^ 0 u 

The. bandwidth of the heading regulator is 
3 rad/sec. 


* 

^20 


^30 


5 

V30 


y3 = E3<';') 


4 

^3 



Y 


5 

3 


The total regulator output, 6y^, is added to Regulator — the heading channel, 

the open-loop command, y^, resulting in the total 

canonical control, y\ it is then transformed u- u ■ ^ 

by means of the W-map given by Eqs. (42) and (43) into the natural control u which, in turn, drives the 

actual plant. 


8. SYSTEM PERFORMANCE 


The results of a manned simulation are summarized in this section. The code was implemented onthe flight 
computer to be used in the flight test, and the mathematical model of the helicopter (UH-IH) was driven 
through the actual hydraulics. 


The experimental flightpath is defined by a set 
of way points, segments of lines and helixes, and a 
speed profile, as shown in Figs. 12 and 13. As can 
be seen from Fig. 12, the flightpath is a closed 
curve. The time dependence is shown in Fig. 13. 

The experiment, which consists of automatically 
flying this trajectory, exercises the systerp over 
a wide range of flight conditions. The helicopter 
is taken from hover (way-point 1 in Fig. 12) to 
high-speed (50 m/sec) accelerating, turning, and 
ascending flight. This input to the system is 
coarse, with a variety of discontinuities. The 
required smoothing is provided by the model servo 
discussed in the preceding section. 

For the data presented, the helicopter was 
flown manually to the point ^ marked in Fig. 12. 
There the automatic system was engaged. It takes 
about 500 sec for the helicopter to go once around 
the flightpath. Unlike the coarse accelerations in 
Fig. 13, yie model accelerations are smooth, as is 
the vertical velocity, Vjq* The second panel in 
Fig. 14 (labeled "acceleration error") shows the 
effects of the neglected parasitic terms on acceler- 
ation. The acceleration errors are quite small - 
less than 0.05 g. The regulator controls these 
effects by means of position errors. The resulting 
horizontal error is less than 2 m, and the vertical 
error is below 0.5 m. The speed error, 

6v ■ II v|i " ll^oll* below 0.5 m/sec. 

Thus, the performance of the design in the sim- 
ulation tests was good, and. at this writing, flight 
tests are in progress. 



Fig. 12. Experimental flightpath shown in 
horizontal and vertical planes. 




ACCELERATION 



Fig. 13. Time-dependence of the coarse command. 



Fig. 14. System response - canonical variables. 
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